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Abstract 



both perturbative QCD and effective field theory in their respective domains of validity and 
that models the transition between them in analogy to the finite temperature QCD phase 



transition. Using techniques of quantum kinetic theory, we formulate a real-time description 
in momentum and position space. The approach is applied to the evolution of fragmenting 
qq and gg jets as the system evolves from the initial 2-jet, via parton multiplication and 
cluster formation, to the final yield of hadrons. We investigate time scale of the transition, 
energy dependence, cluster size and mass distributions, and compare our results for particle 
production and Bose-Einstein correlations with experimental data for e + e~ — * hadrons. 
An interesting possibility to extract the space-time evolution of the system from Bose 
enhancement measurements is suggested. 
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The physics of QCD is well understood in two limits: at short distances (<C 1 fm) 
where the relevant degrees of freedom are quarks and gluons, the dynamics of which is 
accurately described by perturbative QCD, and at large distances ( ^ 1 fm) where the rele- 
vant degrees of freedom are hadrons whose non-perturbative interactions are well described 
by chiral models. What is less understood, and constitutes one of the key open problems 
in QCD, is the transition between the short- and long-distance regimes through interme- 
diate distance scales. This problem is particularly serious for attempts to describe the 
phenomenon of hadronization in high energy QCD processes. Perturbative QCD describes 
quantitatively the short-distance evolution of the dynamical system and hence the large- 
scale flow of energy-momentum. The transformation of colored quarks and gluons into 
colorless hadrons, however, is commonly believed not to involve large momentum trans- 
fers, as expected for a large-distance effect. Moreover, experiments on many high-energy 
QCD processes, including for example e + e~ — ► hadrons and deep inelastic lepton-nucleon 
scattering, strongly support the idea that the fragmentation of partons into hadrons is a 
universal mechanism. However, the theoretical tools currently available for studying QCD 
are inadequate to describe the dynamics of this transformation from partonic to hadronic 
degrees of freedom. Perturbative techniques are limited to the short distance regime where 
confinement is not apparent Q , whilst effective low-energy chiral models JH and QCD sum 
rules Q that incorporate confinement, lack partonic degrees of freedom. In principle, lat- 
tice QCD [Q] should be able to bridge the gap, but in practice dynamical calculations of 
parton-hadron conversion are not yet feasable. 

In this paper we extend previous work Q and advocate a new approach to the hadroniza- 
tion problem that incorporates both perturbative QCD and effective field theory in their 
respective domains of validity and that models the transition between them using ideas 
developed in phenomenological descriptions of the finite temperature transitions from a 
quark-gluon plasma to a hadronic phase ||. The latter is described by an effective theory 
that incorporates a chiral field U whose vacuum expectation value (vev) Uq = (0\U + U^\0) 
represents the quark condensate {0\qq\0), and a scalar field x whose vev xo = (0| x| 0) rep- 
resents the gluon condensate (Ol-F^F^lO). We visualize a high energy collision such as 
e + e~ — > qq as producing a "hot spot" in which the long range order, represented by Uq and 
X0; is disrupted locally by the appearance of a bubble of the naive perturbative vacuum in 
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which (0|gg|0) = = (0\F^ u F^\0). Within this bubble, a parton shower develops in the 
usual perturbative way, with the hot spot expanding and cooling in an irregular stochas- 
tic manner described by QCD transport equations. This perturbative description remains 
appropriate in any phase-space region of the shower where the local energy density is large 
compared with the difference in energy density between the perturbative partonic and the 
non-perturbative hadronic vacua. When this condition is no longer satisfied, a bubble of 
hadronic vacuum may be formed with a probability determined by statistical-mechanical 
considerations. A complete description of this conversion requires a treatment combining 
partonic and hadronic degress of freedom, which is the essential aspect of our approach. 
Although studies of the finite temperature QCD phase transition J?]] indicate that it may 
be completed quite rapidly, in which case mixed description is not needed in a first ap- 
proximate treatment, we are more ambitious here and propose a universal approach to the 
dynamic transition between partons and hadrons based on an effective QCD field theory 
description and relativistic kinetic theory that includes a mixture of both sets of degrees of 
freedom. 

The purpose of this letter is to present the essential concepts of our approach, i.e. the 
quantitative formulation and its application of the above picture. An extensive documen- 
tation of our work can be found in Ref. j|] to which we refer for details. Let us begin by 
defining the distance measure L for the space-time separation between two points r and r' 



and introduce a characteristic length scale L c that separates short distance (L <C L c ) and 
long range (L>L C ) physics in QCD. The scale L c can be associated with the confinement 
length of the order of a hadron radius. In the limit of short distances L <C L c (or high 
momenta, high temperatures), QCD may be described perturbatively by a Lagrangian in 
terms of the elementary gluon (A^ 1 ) and quark fields (tjj,ip), 



where F£ u = A v a — d v A% + g s f a b c A^A" and summation over color indices a, i,j is un- 
derstood. The functions kl and [il introduce an explicit scale(L)-dependence in Cl which 
modifies the quark and gluon properties when L increases towards L c and beyond. In 



[r = r M = (t, r)]: 
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the limit L — > 0, kl = 1 and /u^ = (neglecting the quark current masses), thus the 
QCD Lagrangian is recovered. However, at larger L, the bare quark and gluon fields be- 
come dressed by non-perturbative dynamics and we expect kl < 1 an d [il 7^ 0. In fact, 
the short range behaviour of kl (and similarly of hl) can be calculated perturbatively: 
kl = [1 + 5 f s/(87r) 2 (H — 2nj/3) ln(L A)] -1 . Hence both and the dynamical quark masses 
fiL vary (as does the coupling constant g s ) with the renormalization scale A, which we 
expect to be directly related to the confinement scale L c . We will return to that issue 
below. 

In the limit of large distances L^L C (or low momenta, low temperatures) the hadronic 
physics can well be described by an effective field theory that embodies the scale and chiral 
constraints of the fundamental QCD Lagrangian. The corresponding effective Lagrangian 
is written in terms of collective fields x, U, |q]: 



\{d, X )(d»x) + \tt 



v(x,u) , 



(3) 



where \ is the aforementioned scalar field and U = exp (i J2^=o ^jfy/fir) the pseu- 
doscalar field for the nonet of meson fields (f)j, and the potential V(x, U) is given by 
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which models the quantum dynamics in the low energy regime associated with the self- and 
mutual interactions of x an d U. Here the parameter b is related to the conventional bag 
constant B by B = 6xq/4, c is a constant of mass dimension 3, m q = diag(m M , m^, m s ) is 
the light quark mass matrix. The feu's of the collective fields x an d U are zero in the short 
distance limit L <C L c , i.e. in the naive perturbative vacuum, but become finite in the long 
wavelength limit at large distances L>L C , namely (0|x|0) = xo anci (0|£7 + U^\0) = Uq, 
which are to be regarded as order parameters of the physical vacuum. Thus, the potential 
(H) has a minimum when (x) = Xo an d equals the vacuum pressure (bag constant) B at 

(x) = o. 

As already advertised, we seek to describe parton-hadron conversion in the aftermath of 
a high-energy collision, in which an initial "hot spot" expands irregularly. Inhomogenities 
that appear during the expansion are characterized by local dynamical scales. It is well 



known how to incorporate small distance (high momentum) scales using the renormalization 
group in perturbative QCD. Dynamical scale dependence can also be taken into account in 
the long distance effective Lagrangian (g). For instance, at finite temperature T the effective 
potential (Q) acquires additional T-dependent terms §, resulting in a modification of the 
vev's of x an d U + U> and thus of the gluon and quark condensdates, (Oji^-F^lO) and 
(0|g(7|0). In the present situation we are dealing however with general non-equilibrium 
systems, in which we advocate the length scale L to control the scale dependence of local 
fluctuations. A complete treatment of the hadronization process clearly requires the linkage 
between the elementary partonic degrees of freedom (A^,^) and the collective hadronic 
degrees of freedom (x, U). From the above discussion we expect that 

lim kl = 1 , lim lll = , lim(x)L = 

L^O L^O L->0 

lim k l = , lim /j, L = oo , lim (x}l = Xo (5) 

L^oo L^oo L^oo 

corresponding to the unconfined and confined phases, respectively. In view of the one- 
to-one relation between L and (x)l one may express (||) as boundary conditions of the 
(x)-dependence of kl and when (x) —* or (x) — ► Xo- This characterizes the so-called 
color-dielectric property of the vacuum in QCD which has been widely applied in soliton 
models and hadron phenomenology. 

It is important to realize that the variation of the internal length scale L, as defined 
in (|l[), is governed by the dynamics of the fields itself and in turn it must determine the 
time evolution of the fields A^,ijj as well as x- With regard to (||), we assume that the 
behavior of the functions kl and hl is correlated with the value of (x) and express this 
by kl = kl(x) and hl = Hl{x)- Then, by combining ©-(§), we obtain an effective field 
theory description covering the full range < L < oo, 

C l [A»,tP, X ,U] := C L [A^^\ + £\x,U,tf] , (6) 

in which the elementary gluon and quark fields are coupled to the collective field x via K L (x) 
and hl(x) i n a self-consistent manner that is controlled by the dynamically varying scale L. 
Our present understanding of QCD does not yield us explicit interpolating formulae for kl 
and from first principles. Fortunately, as we find, the particularities of their functional 
forms are not crucial for the treatment of parton hadron conversion that we present here, 



4 



as long as they are smooth and satisfy the boundary conditions (|5|). The reason is that, 
as L decreases, the appropriate dynamical description of QCD apparently jumps rapidly 
from the first limit in (||) to the neighborhood of the second limit, corresponding to the 
rapid formation of hadronic domains when L > L c . This feature is very reminiscent of 
the weakly first order nature of the QCD phase transition at finite temperature. The 
essential point is that, due to the conditions the variation of kl with L and x generates 
color charge confinement at large distances, because of the fact that a color electric charge 
creates a displacement D a = klEo,, where E% = F® k , with energy ^ / d 3 rD 2 / 'kl which 
becomes infinite at large L for non-zero total charge. Similarly, the increase of the dynamical 
mass fj-i with L and X leads to an effective confinement potential that ensures absolute 
confinement also for quarks by prohibiting their propagation at large distances. We tested 
various choices for kl and [Xl and, to be specific, in what follows we adopt the forms 
k l{x) = 1 — (Lx) 2 / (LqXo) 2 , which is a minimal possibility compatible with the QCD scaling 
properties, and Hl(x) = Mo( re £ (x) ~ 1) with constant fXo set equal to 1 GeV. The latter 
form reflects that the quark mass term jil is assumed to be induced by non-perturbative 
gluon interactions (through kl), rather than being an independent quantity, as is suggested 
by an explicit calculation [|l0| of the quark self-energy involving the gluon propagator in 
the presence of the collective field X - 

Let us emphasize once more the similarity of the above approach with finite tempera- 
ture QCD phenomenology ||]. The effect of kl and [il can be interpreted as a scale (L) 
dependent modification 5V, which adds to the (L independent) potential V, eq. @, 

V(L) := V( X ,U) + 5V(L,x) , SV(L tX ) = F^ a F^ + ^(x) , (7) 

with SV = 0(L 2 ). Here the length scale L corresponds to the temperature T in finite 
temperature QCD, where the correction to the zero temperature potential is 0(T 2 ). Hence 
we expect the characteristic confinement length L c to play a similar role as the critical 
temperature T c in the QCD phase transition. However, this formal anology is to be taken 
with some caution, because here we are aiming to describe the evolution of a general non- 
equilibrium system in real time and Minkowski space, as opposed to the thermal evolution 
in Euclidian space. 

We complete our discussion of the above field theory aspects with the following impor- 
tant remarks: (i) The formulation is gauge and Lorentz invariant and is consistent with 
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scale and chiral symmetry properties of QCD. It interpolates from the high momentum 
(short distance) QCD phase with unconfined gluon and quark degrees of freedom and chi- 
ral symmetry ((x) = 0, (U) = 0, Kj_, = 1, = 0), to a low energy (long range) QCD phase 
with confinement and broken chiral symmetry ({x) = %q, (U) = Uq,k,l = 0, hl = oo). 
(ii) By construction (|6|) strictly avoids double counting of degrees of freedom, because the 
introduction of the scale L and the behaviour of the L-dependent the functions kl and \il 
truncate the dynamics of the elementary fields A^, ip to the short distance regime, whereas 
the effective description in terms of the collective fields Xi U covers the complementary long 
range domain, (iii) There is no need for explicit renormalization of the composite fields x 
and U, because they are already interpreted as effective long range degrees of freedom with 
loop corrections implicitely included in V(x, U) and it would be double counting to add 
them again. 

As documented in detail in ||, 11], it is possible to derive from the Lagrangian (g) a 



formulation of QCD transport theory that incorporates both partons and hadrons, yielding 
a fully dynamical description of the QCD matter in real time and complete phase-space. 
Since this serves as the basis for our subsequent description of the conversion of partons 
into hadrons locally in phase-space, we sketch here this formulation of QCD transport 
theory. The first step is to derive from the field equations of motion the corresponding 
Dyson- Schwinger equations for the real-time Greens fuctions of the fields if), A^, x, an d U, 
denoted by G a (x, y), and where a = q, g, x, U. They are defined as the two-point functions 
that measure the time ordered correlations between the fields at space-time points x and 
y. In symbolic notation, one obtains a coupled system of integral equations, 

G a (x,y) = G$\x,y) + W dVdV G^(x, x') E^x', x") G a (x", y) , (8) 

13 J 

where Ga^ denotes the free field Greens functions that satisfy the equations of motion in 
absence of interactions, and the self-energies S Q embody both the mutual and the self- 
interactions of the fields. The explicit expressions can be found in Ref. Q. Next we 
introduce the Wigner transforms 

G a (r,p) = J d 4 Re ip - R G a (x,y) , (9) 

(where r = ^(x + y), and p is the conjugate variable to R = x — y) of the Greens functions 
and similarly for the self-energies S a (x,y). The dependences on r would be trivial for 



field theory in vacuo, but non-trivial in an inhomogenous QCD matter. In terms of the 
Wigner transforms Q a (r,p) and S a (r,p), the Dyson-Schwinger equations (||) become kinetic 
equations describing the real-time evolution of the matter in phase-space spanned by r and 
p^ = (E,p). Then, by tracing over color and spin polarizations, and taking the expectation 
values (or, in medium the ensemble average) of the Wigner transformed Greens functions 
Q a , one obtains the scalar functions 



The c-number functions F a (r,p) for the particle species a are the quantum mechanical 
analogues to the classical phase-space distributions that measure the number of particles at 
time t in a phase-space element d?rcfip. The F a contain the essential microscopic informa- 
tion required for a statistical description of the time-evolution of a many-particle system in 
complete phase-space and provide the basis for calculating macroscopic observables in the 
framework of relativistic kinetic theory. In particular, the local space-time dependent par- 
ticle currents n a for the different particle species and the corresponding energy-momentum 
tensors Tg u are given by [[[^] 



where dVL a = 'y a dM 2 d' i p/(l6iT 3 p ), the 7 Q are degeneracy factors for the internal degrees of 
freedom (color, spin, etc.), M a measures the amount by which a particle a is off mass-shell 
as a result of the selfenergy terms in (||), and p° = E = +y / p 2 + M 2 . These macroscopic 
quantities can be written in Lorentz invariant form by introducing for each species a the 
associated matter flow velocity u^(r), defined as a unit-norm time- like vector at each space- 
time point, {u^u^) a = 1. A natural choice is e.g. -u^ = n£/ 'y 'n v a n v a . By contracting 
the quantities ( p"T| ) with the local flow velocities u^, one can now obtain corresponding 
invariant scalars of particle density, pressure, and energy density, for each particle species 
a individually, 



The above formulation is readily applicable to the dynamics of the parton-hadron con- 
version in rather general situations. Here we will as an illustrative application study the 
fragmentation of a qq jet system with its emitted bremsstrahlung gluons and describe the 




(10) 







(12) 
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evolution of the system as it converts from the parton phase to the hadronic phase. In 
this case the kinetic equations that derive from the Dyson- Schwinger equations (||) simplify 
considerably and yield a set of coupled tranport equations for the phase-space densities 
F a (t,f;p, M 2 ) of the generic form M 

Pli d? + (M a d?M a )di 

processes k 

These kinetic equations reflect a probabilistic interpretation of the evolution in terms of 
successive interaction processes k, in which the change of the particle distributions F a is 
governed by the balance of gain (+) and loss (— ) terms. The left hand left side describes 
propagation of a quantum of species a in the presence of the mean field M a generated by 
the others, and on the right hand side the integral operators incorporate the effects of 
the self-energies in terms of real and virtual interactions that lead to particle excitations. 

We now discuss the illustrative application of the outlined formalism to parton-hadron 
conversion starting from a highly virtual 7* or a Z° in an e + e~ annihilation event with 
large invariant mass Q 2 3> A 2 corresponding to a very small initial Ld fm, that initiates 
a qq jet evolution. At first the partons multiply, determined by a coherent parton shower 
simulation incorporating (coherent) gluon emission, secondary qq production, etc.. Eventu- 
ally, as the quanta diffuse in space-time and L increases, they will coalesce to pre-hadronic 
clusters by tunneling through the potential barrier imposed by the dynamically changing 
V(L), which is determined within a coalescence model. Finally these pre-hadronic clusters 
will convert into physical hadronic states and subsequently decay into low mass hadrons, 
depending on the density of accessible hadronic states. The system of particles is evolved 
in discrete time steps, here taken as At = 0.01 fm, in coarse grained 7-dimensional phase- 
space with cells Afi = A 3 rA 3 pAM 2 . The partons propagate along classical trajectories 
until they interact, i.e. decay (branchings) or recombine (cluster formation). Similarly, the 
formed clusters travel along straight lines until they decay into hadrons. The corresponding 
probabilities and time scales of interactions are sampled stochastically from the relevant 
probability distributions. With this concept, we can trace the space-time evolution of the 
system H: In each time step, any "hot" off-shell parton is allowed to decay into "cooler" 
partons, with a probability determined by its virtuality and life-time. Also in each step, ev- 
ery parton and its nearest spatial neighbor are considered as defining a fictious pre-hadronic 



F„ 



(13) 
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bubble with invariant radius L, as defined by (jlj). Depending on the value of L, the relative 
probabilities for which configuration is more favorable determine whether the partons con- 
tinue in their shower development, or a parton cluster is formed. This cascade evolution is 
followed until all partons have converted, and all clusters have decayed into final hadrons. 

To this end we need to specify the parameters in the potential V, eq. (||), or V(L), eq. 
(0). As L varies, V changes its shape and affects the evolution of the system. The details of 
the dynamics are controlled the choice of bag constant B which defines the vacuum pressure 
V(0) in the short distance limit L — > 0, and xo, the value of the condensation of X lXi the 
long distance regime. Although the values of B and xo are not precisely known, there 
is agreement of various phenomenological determinations about their range: one expects 
HI B 1 / 4 = (150 - 250) MeV and xo = (50 - 200) MeV. In what follows, we choose two 
representative combinations: (B 1 / 4 , X o) = (230, 200) MeV and (B 1 / 4 , X o) = (180, 100) MeV. 
We then proceed in analogy to the QCD phase transition model of Ref. || and compare 
in each time step the local pressure of partons P qg (t,f,L) with the pressure P x (t,r,L) 
that a pre-hadronic bubble would create instead. We compute the pressures from the 
corresponding phase-space densities (|T0|) and the formulae (11), (|T2|) . Representing 



P qg (r,L) = a qg (r,L) L~ 4 — B , P x (r,L) = a x (r,L) L~ 4 - V(L) , (14) 

and, defining L = L c as the characteristic length scale, in analogy to the critical temperature 
in the QCD phase transition, such that the two pressures equal each other at L c and 
V(x,L c ) = 0, we get the determining condition: 

~a qg (r,L c ) - a x (r,L c )^ l! 1 



B 



(15) 



With the functions a qg and a x obtained from the numerical simulation, we find then for the 
characteristic confinement length L c = 0.6 fm for the choice B x l 4 = 230 MeV and L c = 0.8 
fm for B 1 / 4 = 180 MeV. 

Fig. la compares the parton and hadron pressures P(t, L) (calculated at the peak of 
the expanding shock front) as a function of time t in parton showers initiated by (/(/-pairs 
with center-of-mass energies Q = 100 GeV. We see that the pressures cross over at a time 
t ~ 0.7 fm, almost independent of the choice of B (and hence L c ), and we have found 
the same feature in (/(/-initiated showers. The corresponding time development of the total 
generated transverse momentum p±(t) is shown in Fig. lb, where we see that the share 



of hadronic clusters dominates when t £ 1 fm, which again is also true for ^-initiated 
showers. At lower energy (we performed an identical analysis at Q = 10 GeV) we found 
the same features, except that the crossover times are somewhat shorter. 

In this dynamically evolving situation the probability that two partons with invariant 
spatial separation L coalesce to a cluster (viewed as a pre-hadronic bubble formed in the 
vacuum), is determined by ir(L) = 1 — exp(— AF L), where AF is the free energy of a thin- 
walled bubble of one phase (partons), immersed in a medium in the other phase (vacuum): 

AF = — R 2 C a , (16) 

where R c = 2a/ AP is the critical bubble size determined in terms of the surface tension a 
and the difference in pressures AP = P qg — P x . We estimate the surface tension from the 
relation g|] a = a(L c ) = J dx\/2V{L)\ L=Lc with the potential (0), which yields a 1 / 3 = 40 
(48) MeV (for B l / A = 230 (180) MeV). We remark that similar numbers are indicated by lat- 
tice QCD simulations. These small values of a correspond to a weakly first-order transition 
at finite temperature, which is consistent with astrophysical constraints on inhomogenities. 
As a consequence of this feature, although our approach allows for parton-hadron conver- 
sion when the parton separation L is even less than L c , or, when L 3> L c , most clusters 
in fact form when L is only slightly larger than L c , as seen in Fig. 2a. We also observe 
that the distribution of cluster sizes is essentially independent of the initial state energy. 
The same feature was found true for ^-initiated showers. The spectrum of cluster masses, 
determined by the total invariant masses of coalescing partons at the moment of conversion, 
is shown in Fig. 2b which exhibits no sensitivity to the value of L c at all. Moreover, the 
shapes of both spectra turn out to be rather independent of the jet energy, as well as of the 
type of initial partons, and therefore appear to be universal. 

To obtain the hadron spectrom from the cluster distribution, we implement essentially 
the cluster-fragmentation scheme of [|14[, but with the important modification that a cluster 
weighing more than a critical mass M cr u = AGeV is not allowed to convert from partons 
to hadrons, but is forced to continue its partonic shower development. This restriction is 
motivated by the experimental indication that very heavy cluster production is suppressed. 
Although conceptionally significant, it is not very important numerically, because it affects 
less than 5 % of the clusters. We show in Fig. 3a the total charged particle multiplicity n c h, 
which grows with energy in consistency with experimental data. The momentum spectra 
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of charged hadrons with respect to x = 2E/Q, the particle energy normalized to the total 
energy Q, at Q = 34 GeV and Q = 91 GeV, are shown in Fig. 3b and again agree with the 
data. Parton shower Monte Carlos are well known to describe these data sucessfully: our 
repetition of this success simply means that our novel parton-hadron conversion mechanism 
is innocuous in this regard. The comparisons in Fig. 3 do not indicate any clear preference 
for one value of B or L c over the other. 

However, some indication may eventually be drawn from Bose-Einstein correlations 
among produced hadrons. The observed Bose enhancements in the pion spectra are gen- 
erally interpreted as reflecting the size and coherence of the hadron emission region, which 
is directly related to the choice of B and thus L c in our model. We have used the method 
of Sjostrand |18| to simulate the enhancement for same-sign charged pion pairs. We see in 
Fig. 4a that the enhancements are indeed distinguishable for B l l A = 230 MeV (L c = 0.6 
fm) and B 1 ^ = 180 MeV (L c = 0.8 fm). Mor eover, these enhancements are essentially 
independent of Q, as can be seen from the ratios of the enhancement factors plotted in 
Fig. 4b. We compare our predictions with OPAL data in Fig. 4a. Prima facie, these are 
very consistent with our predictions for B 1 ^ = 180 MeV (L c = 0.8 fm). However, the 
unambiguous definition of a preferred value of this input parameter must await a systematic 
analysis of available data on Bose-Einstein correlations in e + e~ — * hadrons. Nevertheless, 
we are pleased that our model appears to be able to correlate sucessfully such diverse quan- 
tities as the bag constant B, the characteristic microscopic scale L c for the parton-hadron 
transition, and measured Bose-Einstein correlations. 

In conclusion: we presented a novel approach to the dynamics of parton-hadron con- 
version and confinement, based on a kinetic multi-particle description in real time and 
complete phase-space, which combines perturbative QCD at high energies and an effective 
field theory approach to hadrons at low energies. As a test application that exhibits generic 
features, we have considered the prototype reaction e + e~ — * hadrons. The main results 
are: (i) the local conversion of partons to hadronic clusters occurs very rapidly, but the 
global time scale for the transition of all parts of the system is comparatively long; (ii) 
features of the perturbative parton evolution are projected unscathed onto hadron distri- 
butions; (iii) the Bose enhancement of same-sign pions may provide a sensitive probe of the 
details of the space-time evolution. Our approach may be extended to other applications 
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involving parton cascades and hadronization in different environments, including deep in- 
elastic lepton-nucleon or -nucleus scattering, high energy hadron-hadron, hadron- nucleus, 
or nucleus-nucleus collisions, and the (non-equilibrium) dynamics of the QCD phase tran- 
sition at finite density and temperature. The development of a systematic understanding 
of the hadron emission regions in different processes appears within reach. 
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FIGURE CAPTIONS 



Figure 1 : a) Time evolution of the kinetic pressures P qg of partons and P x of prehadronic 
clusters for qq initiated jet evolution with total jet energy Q = 100 GeV. b) Total transverse 
momentum p± generated during the time evolution of the system in the center-of-mass of 
the initial two-jet. 

Figure 2: a) Distribution of the cluster sizes of clusters formed from neighboring partons. 
b) Associated cluster mass spectrum. 

Figure 3: a) Resulting average charged multiplicity versus total energy Q in e + e~ anni- 
hilation events, in comparison with experimental data, b) x-spectra of charged hadrons 
(x = 2E/Q) with respect to the variable ln(l/x) at Q = 34 GeV and Q = 91 GeV, 
confronted with measured disributions at PEP and LEP. 

Figure 4: a) Simulated Bose-Einstein enhancement bL c (q) as a function of pair mass q 
of emitted same-sign pion pairs for the two values of L c . b) Ratios of the enhancements 
fy). 6 fm{q)/bo.8 fm(q) f° r total jet energies Q = 34 GeV and Q = 91 GeV. 
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